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Terahertz (THz) lasers on optically pumped multiple-graphene-layer (MGL) structures as their 
active region are proposed and evaluated. The developed device model accounts for the interband 
and intraband transitions in the degenerate electron-hole plasma generated by optical radiation in 
the MGL structure and the losses in the slot or dielectric waveguide. The THz laser gain and the 
conditions of THz lasing are found. It is shown that the lasers under consideration can operate at 
frequencies > 1 THz at room temperatures. 



I. INTRODUCTION 



Graphene, graphene nanoribbons, and graphene bilay- 
ers (see, for instance, Ref. \w) can be used in different 
terahertz (THz) devices @,H 1, H I, 0, II @- Opti- 
cal excitation of graphene can result in the interband 
population inversion and negative real part of the dy- 
namic conductivity of a graphene layer in the THz 
range of frequencies @, H, H|. The negativity of the 
real part of the dynamic conductivity implies that the 
interband emission of photons with the energy fku, where 
H is the reduced Planck constant, prevail over the inra- 
band (Drude) absorption. If the THz photon losses in 
the resonant cavity are sufficiently small, the THz lasing 
can be realized in graphene-based devices with optical 
pumping [l(| EH • In Refs. [l(| E| , the optically pumped 
THz lasers with a Fabri-Perot resonator were considered. 
In this paper, we propose and evaluate the THz lasers 
based on multiple-graphene-layer (MGL) structures with 
a metal slot-line waveguide (SLW) or a dielectric waveg- 
uide (DW) pumped by optical radiation. The specific 
features of characteristics of the MGL-based lasers un- 
der consideration are associated with the frequency de- 
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pendences of the absorption in the waveguides and the 
gain-overlap factor, which is sensitive to the spatial dis- 
tribution of the THz electric field. As show in the follow- 
ing, the characteristics of the lasers with SLW and DW 
are fairly similar, so that we will mainly focus on the de- 
vice with SLW. This in part because the structures with 
SLW exhibit somewhat better confinement and can be 
used not only for MGL-based lasers with optical pump- 
ing but also in MGL-based injection lasers with lateral 
p-i-n junctions (both electrically induced [H,[ll,[l3| and 
formed by pertinent doping of MGL structure). 



II. DEVICE MODEL 

We consider lasers with an MGL structure on a SiC or 
Si substrate with the side highly conducting metal strips 
and a highly conducting back electrode at the substrate 
bottom or on the top of a DW. The cross-sections (cor- 
responding to the y — z plane), of the device structures 
under consideration are schematically shown in Fig. 1. 
The axis x correponds to the direction of the electromag- 
netic wave propagation, whereas the y and z directions 
are in the MGL structure plane and perpendicular to it, 
respectively (see Fig. 1). The MGL plane corresponds 
to z = 0. The finiteness of the MGL structure thickness 
can be disregarded. It is assumed that the MGL struc- 
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FIG. 1: Schematic views of the device structures under con- 
sideration: (a) with a SLW and (b) with a DW. 
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ture under consideration comprises K upper GLs and a 
highly conducting bottom GL on a SiC substrate or K 
GLs (without the bottom GL) on a Si substrate. Epi- 
taxial MGL structures with up to K = 100 GLs with 
very long momentum relaxation time of electrons and 
holes (r ~ 20 ps) were recently fabricated using the ther- 
mal decomposition from 4H-SiC substrate [la]. MGL 
structures without the bottom GL can be fabricated us- 
ing chemical/mechanical reactions and transferred sub- 
strate techniques which include chemically etching the 
substrate and the highly conducting bottom GL [Ty] (or 
mechanically peeling the upper GLs) and transferring the 
upper portion of the MGL structure on a Si or equiva- 
lent transparent substrate. Since the electron density 
and Fermi energy e% in the bottom GL is rather large 
(ef ~ 400 meV [Hi), the Drude absorption in this GL 
can be significant although it can be overcome by a strong 
emission from the upper GLs if their number is suffi- 
ciently large. The MGL structures without the bottom 
GL can exhibit significant advantages (see Ref. [ll[ and 
below) . 

It is assumed that the MGL structure is illuminated 
from the top by light with the energy of photons htt. 
The optical waveguide input of the pumping radiation 
is also possible. When M7 is close to Nhuio/2, where 
hojQ ~ 0.2 eV is the optical phonon energy and N is an 
integer, the photogeneration of electrons and holes and 
their cooling associated with the cascade emission of op- 
tical phonons, result in an essentially occupation (popu- 
lation inversion) of low energy states near the bottom of 



FIG. 2: Coefficients of (a) absorption and (b) reflection versus 
frequency for SLWs with different geometrical parameters. 



the conduction band and the top of the valence band. At 
elevated electron and hole densities (i.e., at sufficiently 
strong optical pumping), the electron and hole distribu- 
tions in the range of energies e <C hiog in the fc— th GL 
(1 < k < K) can be described by the Fermi functions 
with the quasi- Fermi energies . 

The quasi-Fermi energies in the GLs with k > 1 are 
mainly determined by the electron (hole) density in this 

layer i.e, ey oc and, therefore, by the rate 

of photogeneration G n by the optical radiation at the 
k— th GL plane. Considering the attenuation of the opti- 
cal pumping radiation due to its absorption in each GL, 
one can obtain 



= - l3) K - k + (1 - Pb?(1 [3) K+k - 1 ]. (1) 

Here Iq is the intensity of incident pumping radiation, 
(3 = ne 2 /he ~ 0.023, where e is the electron charge, c 
is the speed of light in vacuum, and 0g — (4tt / c)Kea^ . 
The latter quantity accounts for the absorption of optical 
pumping radiation in the bottom layer. A relationship 
between and Gq ^ is determined by the recombina- 
tion mechanisms As a result, £p can be expressed 
via the quasi-Fermi energy in the topmost GL ep" = e^P , 



3 




IEJ, arb. units 





0.8750 
1 1.750 

I 2.625 
3.500 
4.375 
5.250 
■ 6.125 
1 7.000 



-20 20 
Coordinate y, urn 




° 

0.8750 
1 1.750 
2.625 
3.500 
4.375 
5.250 
I 6.125 
1 7.000 



-20 20 
Coordinate y, |jm 



FIG. 3: Spatial distributions of THz electric field (a) in SLW 
cross-section and (b) in DW cross-section. 



which, in turn, is a function of the the intensity of inci- 
dent pumping radiation Iq. 

Since the thickness of the MGL structure is small in 
comparison with the wavelength of THz radiation, the 
generation and absorption are determined by the real 
part of the net dynamic conductivity which is the sum of 
the real parts of the dynamic conductivity of the bottom 
GL Re a , and other GLs er^ : 
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Re a u = Re + Re cr 



(2) 



fc=i 



Considering the expressions for Re and Reci/ ob- 
tained previously [111 ], one can arrive at the following: 
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FIG. 4: Frequency dependences of THz gain in laser struc- 
tures with different width of the slot D and different number 
of GL K at different values of the quasi-Fermi energy ej: (a) 
D — 10 /im and (b) D — 15 (jm. 



Here tb and r are the electron and hole momentum relax- 
ation times in the bottom and other GLs, respectively, 
T is the electron and hole temperature, and ks is the 
Boltzmann constant. The first two terms in the right- 
hand side of Eq. (3) are associated with the intraband 
(Drudc) absorption of THz radiation in all GLs, whereas 
the third term corresponds to the interband transitions. 
In the case of lasers without the bottom GL, one can use 
Eq. (3) formally setting tb — oo. 
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FIG. 5: The same as in Fig. 4 but for laser structures without 
bottom GL. 



The electromagnetic waves propagating along the SLW 
was considered using the following equation which is a 
consequence of the Maxwell equations: 



d 2 EUy,z) , d 2 E L0 { Vl z) 



r w 2 



dy 2 



dz 2 



E u {y,z) = Q. 

(4) 

Here E^y, z) the amplitude of the y— th com- 
ponent of the THz electric field E(t, x, y,z,) — 
E ul (y, z) exp[i(qx — wt)], rj(y,z) is the complex per- 
mittivity , q is the wave number of the propagat- 



ing mode. The quantities E w (y,z), dE UJ (y, z) /dz, and 
ij~ 1 (y, z)d[rj(y, z)E UJ (y, z)]/dy are continuous at the in- 
terfaces between the layers with different refractive in- 
dices. The boundary conditions for the guided mode 
correspond to the condition E u (y, z) — ► at y, z — ► ±oo. 
Equation (4) was solved numerically using the effective 
index and transfer-matrix methods (see, for instance, 
Refs. [Hill). The coefficient of absorption of the propa- 
gating mode ctuj and the coefficient of reflection from the 
interfaces between the laser structure edges and vacuum 
R were calculated using the following formulas, respec- 
tively: = 21m q and R = \(q — qo)/(q + qo)\ 2 , where 
qo = u Jc. The THz gain, which describes the attenuation 
or amplification of the propagating mode, under optical 
pumping can be calculated using the following formula: 



where 



J%|£Uy,0)| 2 rfy 

I^IZolEUv^Wdydz 



(5) 



(6) 



is the gain-overlap factor (in the case when a u is inde- 
pendent of coordinate y) and r]$ is the permittivity of 
the substate (SiC or Si). 



III. RESULTS AND DISCUSSION 



Figure 2 shows the absorption coefficient a. 



omega 



Of 

electromagnetic waves propagating along the SLW and 
the coefficient of reflection from the laser structure edges 
as function of frequency lo/2tt calculated using Eq. (4) 
for different widths of the slot D and the strips D c and 
the substrate thickness W . The strips and the back elec- 
trode are assumed to be made of Al. Examples of the 
spatial distributions of the THz electric field in the elec- 
tromagnetic wave in a laser structures with D = 15 /im at 
oj/2tt = 1.8 THz in question are show in Fig. 3. The elec- 
tric field distributions with the frequency dependences 
of the dynamic conductivity of the MGL structures ob- 
tained by solving Eq. (4) were substituted to Eqs. (5) 
and (6) to find the gain-overlap factor and the THz 
gain. Figures 4 and 5 show the frequency dependences 
of the THz gain calculated for the laser structures 
with (Fig. 4) and without (Fig. 5) bottom GL and with 
different structural parameters at different pumping con- 
ditions (different values of the quasi-Fermi energy) and 
at T = 300 K. ft was assumed that htt = 920 meV, 
Tg = 1 ps, and r = 1 ps. As seen from comparison of 
Figs. 4 and 5, the laser structures with and without the 
bottom GL exhibit qualitatively similar frequency depen- 
dences of the THz gain g u (at chosen geometrical param- 
eters) but with somewhat higher maxima in the latter 
structures. The frequency u> m in a t which changes its 
sign is uj m i n /2TT ~ 2 THz in the laser structures of both 
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FIG. 6: Comparison of THz gain versus frequency depen- 
dences in laser structures(a) with and (b) without bottom 
GL and with different spacing between MGL structure and 
back electrode W. 
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FIG. 7: The same as in Fig. 6(a) but for K = 100 and e T F 
50 meV. 



types. The latter is different from the situation in the 
MGL lasers with the Fabri-Perot resonator [ll[ in which 
eliminating the bottom GL can lead to a pronounced de- 
crease in uj min . This is attributed to the effect of SLW: 
the variation of the frequency results in a redistribution 
of the spatial distribution of the THz electric field and, 
hence, in a change in the gain-overlap factor. Due to 



FIG. 8: Frequency dependences of THz gain in laser struc- 
tures (without the bottom GL) with different r. 



this, the spacing between the MGL structure W and the 
back electrode affects the the frequency dependences of 
the THz gain. Figure 6 shows these dependences cal- 
culated for different W. One can see that the height 
of the THz gain maxima decreases with increasing W. 
Simultaneously the frequency oj m i n shifts toward lower 
values. This shift is more pronounced in the laser struc- 
tures without the bottom GL [compare Figs. 6(a) and 
6(b)]. As seen from Fig. 6(b), w m i„/27r can reach 1 THz. 
In the case of larger K and , as seen from Fig. 7, the 
THz gain in the range ui/2ir > 1 — 2 THz can be rather 
large. 

A decrease in the momentum relaxation time t (which 
results in an enhancement of the Drude absorption) leads 
to an increase in uj m i n and to a decrease of |Rerj w | in the 
range of frequencies where Re < 0. This is seen in 
Fig. 8. 

The above results correspond to T = 300 K. Lower- 
ing of the temperature should lead to widening of the 
frequency range where Re is negative and where g u 
can be positive. As a result, to m i n might decrease with 
decreasing temperature. The latter is confirmed by the 
calculated temperature dependences of Lo m in shown in 
Fig. 9. 

Figure 10 shows the frequency dependences of the THz 
gain calculated for laser structures with SLW and DW for 
T = 300 K and r = 10 ps. One can see that the maxi- 
mum of the THz gain in the laser with SLW is somewhat 
higher than that in the laser with DW. This can primar- 
ily be explained by the effect of the spatial distribution 
of the THz electric field on the gain-overlap factor [com- 
pare Figs. 3(a) and 3(b)]. To maximize the THz gain at 
the desirable frequency one needs carefully optimize the 
geometrical parameters. 

Considering the propagation of the THz electromag- 
netic wave in the laser structure (with the length L in 
the x-direction) and its reflection from edges or from the 
external mirrors (with the reflective coefficients R± and 
R2), the condition of lasing can be presented as 
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(7) 



For a laser structure with SLW with D — lb fim, W = 
12 /xm, and K = 100 setting ej- — 30 — 50 meV and 
w/2tt = 3.0 THz, so that R x = R 2 ~ 0.05 and g w ~ 
30 — 60 cm -1 (see Figs. 2(b) and 4(b), respectively), one 
obtains L > 0.1 - 0.2 cm. If £> = 15 /xm, IF = 25 /xm, 
and if = 100 at e T F = 50 meV and lu/2tt = 1.5 THz 
(i?i = i?2 — 0.05 and ~ 10 cm" 1 as follows from 
Figs. 2(b) and 7, respectively), one obtains L > 0.6 cm. 

As shown above, at the quasi-Fermi in the topmost 
GL about Ep = 30 — 50 meV, the achievement of the 
THz lasing in the devices under consideration at room 
temperatures is feasible. At T = 300 K the condition 
£ ? ^ 30 meV corresponds to the electron and hole densi- 
ties about 2 x 10 11 cm -2 . Such densities can be obtained 



at reasonable optical powers (see Ref. 
ences therein). 
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FIG. 9: Temperature dependences of uimin for laser structures 
with different r: (a) W = 12 fj,m and (b) W = 25 /xm. 



IV. CONCLUSIONS 
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We proposed THz lasers with optical pumping based 
on MGL structures with waveguides and calculated their 
characteristics. The feasibility of lasing in the devices 
under consideration at the low end of the THz frequency 
range at room temperatures was demonstrated. 
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